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ABSTRACT 

Many exoplanets have now been detected in orbits with ultra-short periods, very close to the Roche 
limit. Building upon our previous work, we study the possibility that mass loss through Roche lobe 
overflow (RLO) may affect the evolution of these planets, and could possibly transform a hot Jupiter 
into a lower-mass planet (hot Neptune or super-Earth). We focus here on systems in which the mass 
loss occurs slowly (“stable mass transfer” in the language of binary star evolution) and we compute 
their evolution in detail with the binary evolution code MESA. We include the effects of tides, RLO, 
irradiation and photo-evaporation of the planet, as well as the stellar wind and magnetic braking. 
Our calculations all start with a hot Jupiter close to its Roche limit, in orbit around a sun-like star. 
The initial orbital decay and onset of RLO are driven by tidal dissipation in the star. We confirm 
that such a system can indeed evolve to produce lower-mass planets in orbits of a few days. The RLO 
phase eventually ends and, depending on the details of the mass transfer and on the planetary core 
mass, the orbital period can remain around a few days for several Gyr. The remnant planets have a 
rocky core and some amount of envelope material, which is slowly removed via photo-evaporation at 
nearly constant orbital period; these have properties resembling many of the observed super-Earths 
and sub-Neptunes. For these remnant planets we also predict an anti-correlation between mass and 
orbital period; very low-mass planets ( M p \ < 5 M®) in ultra-short periods (P or b < 1 d) cannot be 
produced through this type of evolution. 

Subject headings: Planetary Systems: planet-star interactions-planets and satellites: gaseous planets- 
stars: evolution-stars: general-(stars:) planetary systems 


1. INTRODUCTION 

Hot Jupiters, giant planets in orbits of a few days, 
constitute one of the many surprises of exoplanet 
searches. Whether their tight orbits are the re- 
sult of inward migra ti on in a protop l aneta r y disk 
dGoldreich fc Tremainel 119801 : iLin et all 1 19961 iWardl 
119971 iMurrav ct all 1 1 998tb ~ or tidal circularization of 
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Naeasawa et al.1 20081 1 Wu fe Lithwickl 20111 iNaoz etHl 


201H 


Plavchan fe Bilinskil 120131: lYalsecchi fe Rasic 


2014al rbr> . is still matter of debate. Certainly, indepen¬ 


dent of the formation mechanism, tidal dissipation in 
the slowly-spinning host stars is causing the orbits of the 
tighte st hot Jupiters cur r ently known to s hrink rapidly 
(e.g., iRasio et al.l 11996 : iSasselovI 120031 : iBirkbv et ail 
2014 lYalsecchi fe Rasic 1201-lal see also Table 1 in 
Valsecchi fe Rasioll2014a and references therein). 

Eventually, hot Jupiters may decay down to their 
Roche-limit separation. While it is commonly as¬ 
sumed that the planet is then quickly accreted 
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by the star (e.g.. iJackson et al.1 120091: iMetzger et al.1 


2012tlSchlaufman fc Winnll20 1 4 Damiani fc Lanzall20l4 


Teitler fc Koniglir20l4 IZhang fc Pcnevll2014lh for a typ¬ 

ical system hosting a hot Jupiter orbiting a sun-like 
star the ensuing mass transfer (hereafter ‘MT’) may 
be dynamically stable (iScpinskv et a.1.1 [2(TlT) J. This was 
suggested, e.g., t o explain WA S P-12’ s transit features 
(lLai et al.l 120101 '). iTrilling et al.l (jl998H investigated the 
possibility of halting inward disk migration through tides 
and Roche-lobe overflow (‘RLO’) MT from a hot Jupiter 
to a rapidly-spinning (young) stellar host. However, we 
note that the host stars of the tightest hot Jupiters are 
obse rved to be rotat i ng slo wly at present. 

In IValsccchi et al.1 (l20l4h we investigated the fate of 
a hot Jupiter transferring mass to its stellar host us¬ 
ing a simplified binary MT model. We showed that 
the planet could be stripped of its envelope, result¬ 
ing in a hot super-Earth-type planet. This model 
naturally solves some of Kepler ’s evolutionary puzzles 
(e.g., Kepler-78; iHoward et al.1120131 iPene et ahl 120131 
iSanchis-Oieda et al.l 12013T) . and it could explain the ex¬ 
cess of isolated hot super-Eart h- or sub-Neptune-s ize 
planets seen in the Kepler data (jSteffen fc Farrll2013 f). 

Our previous work relied on several simplifying and 
potentially key assumptions. First, while using detailed 
models for the host star, we adopted published mass- 
radius relations for the planetary component, thus as¬ 
suming that the planet remains in thermal equilibrium 
throughout the MT. Second, even though our plane¬ 
tary models included the effect of stellar irradiation on 
the planetary mass-radius relations, irradiation was kept 
fixed, while it is expected to vary as the orbit evolves dur¬ 
ing MT. Furthermore, we neglected the resulting mass 
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loss due to photo-evaporation (IMurrav-Clav et alJl200il 
Jackson et al.l 120101: ILopez et alJl2012t iLopez fc Fortnevl 
20131 lOwen fe Alvaredl2015D . even though various stud¬ 
ies found it to play an important role in the evolu¬ 
tion of highly irradi ated super-Earth and sub-Neptune- 
type planets (e.g. , Uackson et al.1 120101: iRogers et, all 
20111: ILopez et all 120121 iBatvgin fe Stevensonl 120131 
Lopez fe Fortnevll2013l: IQ wen fc Wull2013H~ 

As a natural continuation of our previous work, 
here we significantly expand upon our simple binary 
MT model to include detailed planetary evolution, 
as well as the effects of a varying irradiation and 
the consequent photo-evaporative mass loss from the 
planet. For these new calculations we utilize the Mod¬ 
ules for Experi ments in Stellar As t rophysics ( MESA) 
evolution code (iPaxton et all 1201 lL 120131 201 . The 

MESA inlist files used in this work can be downloaded at 
https: //github.com/FrancescaV/planet _MT_with_MESA 
_inputs_and_models. 

The plan of the paper is as follows. In Section §[2]we 
describe the stellar and planetary models used in this 
work. In Section § [3] we describe our orbital evolution 
model and the physical mechanisms entering the calcu¬ 
lation. We present some examples of our orbital evolu¬ 
tion calculations in Section [I] and discuss our results in 
Section [5] We conclude in Section [G] 

For quick reference, the notations adopted in this work 
are summarized in Table [I] 


Table 1 

Definition of the main parameters used in this work. 


Parameter 

Definition 

M,, M p i 

Mass 

M c , -Menv 

Planetary core and envelope mass 

i?*, R p \ 

Radius 

-^lobe,*? -^lobe,pl 

Roche-lobe radius 

T*, T p i 

Surface temperature 

, ^pl (f^o) 

Spin (orbital) frequency 

z 

Metallicity 

a 

Semimajor axis 

-^orb 

Orbital period 

^orb 

Orbital angular momentum 

f — Menv/Mpi 

Envelope mass fraction 

q = M pl /M* 

Planet to star mass ratio 

t, t MS 

Stellar age and main-sequence lifetime 


Note. — The subscripts and “pi” refer to the star and 
planet, respectively. We take the stellar age to be equal to the 
system age, and the main-sequence lifetime for solar mass stars 
with radiative cores to be the age when the mass fraction of H at 
the center of the star —> 0. 

2. STELLAR AND PLANETARY MODELS 

The stellar and planetary models adopte d in this work 
are computed wi th MESA (version 7184; IPaxton et al.l 
I2011L120131120151) . In particular, the planets are created 
and evolved closely following the test suite make-planets 
provided within MESA and t he input files yielding Fig¬ 
ure 3 of IPaxton et al.l (12013T FI. In what follows we give 

5 available at http : / jmesastar.org / results/mesa2/planets 


specifications only for those MESA parameters that are 
changed from the values adopted in these input files. 

In all our calculations, we consider a 1M 0 star 
paired with a 1 Mj planet. We assume solar com¬ 
position (Y = 0.27, Z = 0.02) for both the star 
and planet envelope, and keep the mixing length pa¬ 
rameter «mlt to MESA’s default value of 2. For the 
plane t, we expand on our previous work (IValsecchi et al.l 
I2014H and consider models with solid cores of masses 
M c = 1 M®, 5 M®, 10 M®, 15 M® and, 30 M®. For the 
cores, we use a constant de nsity of 5gcm -3 , following 
IBatvein fc Stevensonl (|2013D . Furthermore, we take the 
heat-flux arisin g from radioac t ive de cay in the core to be 
zero, following IFortnev et al.l (120071) . As these authors 
point out, this is a common assumption in evolution¬ 
ary models of Jupiter- and Saturn-type planets ( Hubbardl 
i 19771 : fSaumon et alJl!992b IFortnev fc Hubbardn2003f h as 
this approach introduces a small error compared to other 
unknowns entering the problem. However, we note that 
I Lonez fc Fortnevl (120141) found such heating to play an 
important role in delaying cooling and contraction, par¬ 
ticularly for planets less than 5M®. This could lead to 
an underestimate of the radii of sub-Neptune planets, 
especially at ages < 1 Gyr. Below we focus on planetary 
models at least 2 Gyr old. 

We account for the effects of irradiation and photo¬ 
evaporation as follows. For irradiation, we use the E* — 
E p i surface heating function. Here E* is the day-side flux 
from the star at the substellar point, given by 



The planet equilibrium temperature, T eq , is 

/» \ X / 2 

= T * (iJ > ^ 

dSaumon et al.l I1996D so that the power received from 
the host star could be radiated away in equilibrium if 
the planet had this temperature over its entire surface. 
The parameter E p i is the column depth reached by irra¬ 
diation. Here we adopt a value of E p i = lg/cm 2 , which 
yields planetary ma ss-radius relati ons in agreement with 
detailed models by IFortnev et al.l (120071) . within a few 
percent for highly irradiated planets, as shown in Fig¬ 
ure [p W e discuss the effect of varying E p i in Sec¬ 
tion 14.1.31 In what follows, the planetary radius cor¬ 
responds to an optical depth r = 2/3. 

Irradiation leads to photo-evaporative (‘PE’) mass loss 
from the planet. Specifically, this process is thought to 
be most efficient when a hot Jupiter is strongly irra¬ 
diated by ultraviolet (UV) and X-ray photons, which 
photo-ionize atomic H in the planetary atmosphere. 
The resulting heat input, when high enough to induce 
temperatures c orresponding to the escap e velocity, can 
cause outflows. IMurrav-Clav et afl (I2009D identified two 
regimes, based on the stellar flux Fxuv (however, see 
lOwen fc Alvarezll2015t ). For large Fxuv, like those typi¬ 
cal of T-Tauri stars (> 10 5 erg cm -2 s -1 ) the mass loss 
is “radiation/recombination” limited and it is described 
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Figure 1. Planetary mass-radius relations at 4.5 Gyr for different 
core masses and equilibrium temperatures. From top to bottom: 
core-less planets, M c = 10 M®, and 2 5 M®, as illustrative exam¬ 
ples. The solid lines are the models bv I Fortney et ahl <120071 ) (from 
Table 4 of their paper) at an equilibrium temperature of 78 K (light 
grey), 1300K (dark grey), and 1960K (black). The data points 
are the models computed with MESA for E p i set to 1 g/cm 2 (“x") 
and 100g/cm 2 (“□”). Varying E p i between 0.1—1 g/cm 2 makes 
no significant difference. The value of E p i is not im portant for 
T eq < 1300A (symbols overlap). As in I Fortney et all (120071 ). the 
radii correspond to a pressure of bar. We note that the initial 
value of the planetary rad ius required by the createjinitial.model 
routine (see Section 2.1 of I Paxton et al .1120131) was varied between 
2—5 Rj to facilitate MESA’s convergence. 


by 


Mrr-lim ~ 4 X IQ 12 


Fx uv 


5 x 10 5 erg cm -2 s~ 


1/2 


gs 


(3) 

For lower Axuv values the mass loss i s “en e rgy limited” 
and it is described by (lErkaev et al.l 120071 : iLopez et al.1 


120131 


Me -lim 


e7rExuv-Rxuy 
G A I p A tide 


(4) 


where Ar tide = 1 - (3/2)(l/£) + (l/2)(l/£ 3 ) and 

£ = Arm/A xuv- Since in all our calculations 

Ax uv (computed as described below) remains below 
10 5 ergcm~ 2 s _1 , we use only Eqn. 0- However, we 
discuss whether our results are se nsitive to the photo¬ 
evaporation prescription i n Section 14.1.41 an d Tabic [2] 
For the flux, we follow iRibas et al.l ( 20051) and take 


( , \ —1.23 _2 

Gyr) «e s " lc “' 2 ' < 5 > 

where we have scaled their result at 1 AU to an arbi¬ 
trary distance. The parameter e represents the efficiency 
of converting Axuv into usable work, while Axuv is the 
radius of the planet at which the a tmosphere becomes 
optic ally thick to XUV photons. iMurrav-Clav et ahl 
(2009) place Axuv at a surface pressure of ~ 10~ 9 bar, 
which corresponds to a radius that is ty pically 10%—20% 
great er than the optical photosphere (jLopez fc Fortney! 
l20ll hereafter LF13). Here we closely follow the re¬ 
cent results of LF13 and adopt e = 0.1. We note, 
however, that hydrodynamic calculations suggest that e 
can vary b etween 0. 0 1—0.2 depen ding on the mass of 
the planet (lOwen fe Jacksonl 120121) . With e = 0.1, the 
Axuv value that better matches the results of detailed 
numerical calculations by LF13 for systems 1—10 Gyr 
old is Axuv = 1-2 A p i (see below) and we use this es¬ 
timate of Axuv throughout our calculation^. Finally, 
Xtide is t° account for the fact that the mass leaving 
the planet needs only to reach the Hill radius to escape, 

where Anm = M^( 3 (3M*) -1 / 3 a (lErkaev et al.ll2007lh 
We test the implementation of photo-evaporation by 
comparing with the mass loss calculations presented by 
LF13 for the planet Kepler-36c. These are shown in 
Figure [2j Following their work, we create a 9.4 M®, 
lOMyr old irradiated planet with a core of 7.4 Mm a nd 
Z = 0.35. We adopt T eq = 930 K (|Carter et al.l 12012 re- 
ports T eq = 928 ± 10 K) and evolve the planetary model 
with irradiation (fixed) and photo-evaporation, accord¬ 
ing to Eqn. 0- The parameters adopted in this work are 
those yielding agreement with LF13 within a few percent 
for 1 — 10Gyr old systems (black solid line). 


3. ORBITAL EVOLUTION MODEL 

MESA allows us to track the evolution of both the 
planet and star simultaneously, while orbital evolution 
is followed by taking into account changes to the or¬ 
bital angular momentum of the system. Below we de¬ 
scribe the physical mechanisms included in our model: 
photo-evaporatio n (‘PE’), tides, magnetic braking (‘MB’; 
lSkumanichlll972l ). and RLO. For simplicity, we neglect 
the effects of stellar wind mass loss. From the wind pre¬ 
scription provided in MESA’s tes t suite lM_yre-ms-to-wd 
for the evolution of a 1 Mq star (Blocckcr 19951: iReimersI 

6 MESA uses automatic mesh refinement, thus adjusting the 
number of mesh points of the planetary model at the beginning 
of each timestep, if necessary. For this reason, the point at the 
surface where the pressure is ~ 10 <:i bar is not always resolved. 
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Figure 2. Mass loss evolution of Kepler-36c. Percent of mass 
contained in the envelope as a function of time for different Rxuv 
and e values. The grey data points are taken from Figure 1 of LF13. 
The parameter S p i was set to lgcm -2 . For Fxuv? we follow 
I Ribas et al.l (120051) , but rescale the flux at the Kepler-36c orbital 
separation (a = 0.12 AU). As in LF13, we keep -Fxuv constant at 
the 100 Myr value when the star is younger than 100 Myr and we 
let it evolve when the star is older than 100 Myr. Here we assume 
a 1 Mq companion. At the surface (r = 2/3), the pressure is 
lOmbar (where LF13 place the transiting radius). 


fl97l . we find that the orbital evolution timescales asso¬ 
ciated with stellar winds are longer than 10 12 yr through¬ 
out the main-sequence evolution of the host star. 


3.1. Photo-Evaporation 

Mass escape via photo-evaporation affects the orbital 
separation and the spin of the planet. For simplicity, we 
assume spherically symmetric mass loss, which carries 
away the specific angular momentum of the mass losing 
component. With this assumption, photo-evaporation 
affects the orbital angular momentum J or b according to 


( ■/orb \ _ -Mpi'PE 1 _ -Mpi'PE , . 

■/orb J pE -Mpi 1 + q A/pi 

where q = For reference, the evolution of the 

orbital separation due to photo-evaporation is given by 


_ -Mpi.PE q 
\aJ PE M pl 1 + q 


(7) 


Photo-evaporation is always active and Mp^pE is given 
by either Eqn. m or Eqn. depending on Fxuv 
though, in practice, only Eqn. [4] is needed in our calcula¬ 
tions. As far as the planetary spin is concerned, photo¬ 
evaporation carries away the angular momentum of the 
corresponding shells of material. 

While we assume that mass loss is spherically sym¬ 
metric, we note that strong magnetic fields may con¬ 
fine th e flow primarily to the poles and day-side of the 
planet (jOwen fc AdarmJfeOldH . iTbvssandier et al.l (|2015f ) 
studied the torque on super-Earth and sub-Neptune-type 
planets due to anisotropic photo-evaporative mass loss 
using steady-state one-dimensional wind models. They 
found that only in rare cases is the planet’s orbit affected 
by wind torques. 


3.2. Tides 

Tides affect the stellar and planetary spins, as well as 
the orbital separation, transferring angular momentum 
between the orbit and the components’ spin. Within 
MESA, tides are first applied to the spins. The orbital 
separation is then varied so as to conserve total angu¬ 
lar momentum. Here we take the co mponents to be 
rotating as solid b odies (however, see iStevensonl 119791 
IBarker et, al.ll2014ll . 

F or stellar tides we proceed a s in lValsecchi et all (|2014fl 
and IValsecchi fc Rasiol (l2014al fbt) . Specifical ly, we adopt 
the weak-friction approximation (lZahnlll977L ll989D using 
a parametrization for tidal dissipa tion calibrated from 
observations of stellar binaries, as in lHurlev ct alJ ( 2002 1. 
This assumes that tides are dissipated in the stellar con¬ 
vection zone via eddy viscosity. Furthermore, we reduce 
the efficiency of tides at high tidal forcing frequencies 
(when the forcing frequency is higher than the convective 
tur nover f requen cy of the largest eddies) linearly, follow¬ 
ing iZahnl (11966D, and as sug gested by recent numerical 
results bv lPenev et al. ( 200 7). 

For planetary tides, we assume they efficiently main¬ 
tain the planet in a tidally locked configuration (fl p i = 
fl 0 ) throughout the evolution. In fact, for a tidal quality 
factor Q' = 10 6 (typical for gas giants) and a nearly syn¬ 
chronized plane® the spin synchronization timescale due 
to static tides in the planet is ~ 2—4 orders of magnitude 
shorter than the timescale related to the main driver of 
the orbital evolution (i.e., mass loss from the planet; see 
Section ®, depending on the core mass. Clearly, tides 
would synchronize the planet even faster for lower values 
of Q' more appropriate for rocky planets. We further 
discuss tidal locking for the planet in Section [5] 

Even though our calculations account for tides in both 
components, our results show that only stellar tides can 
affect the orbital separation significantly. In particular, 
for a slowly spinning stellar host, tides transfer angular 
momentum from the orbit to the stellar spin, causing 
orbital decay. 


3.3. Magnetic Braking 

For the loss of stellar spi n angular mome ntum via mag¬ 
netic braking we follow iSkumanichl ( 1972 1 and adopt 

(A)mb = (8) 

wher e umb = 1-5 x 10~ 14 yr le.g.. IBarker fc Qgilvid 
20091: iDobbs-Dixon et all 120041 iMatsumura et al.l 120101 : 
Valsecchi fe Rasioll2014al lbh. This law is well established 
for the stellar equatorial rotation rates of interest here 
(1-30 km s” 1 ). 


3.4. Roche-Lobe Overflow 

RLO is modeled by implicitly computing the mass 
transfer rate that is required for the planetary radius 
to remain below its Roche lobe radiu s, for which we use 
the approximation bv lEggletonl (11983H. This proce dure is 
described in Section 2.3.2 of IPaxton et al. i (12015( 1. Here 
we consider both conservative and non-conservative MT. 

7 To get a sense for the magnitude of the s pin synchronization 
timescale we use Eqn. (10) in IMatsumura et all <2010 1 and f2 p i(t) = 
f 1 a (t — At), where A t is the time interval between two consecutive 
time steps during an orbital evolution calculation. 
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When MT is conservative, all mass leaving the planet via 
RLO is accreted onto the star and M* = — M p i ,rlo- In¬ 
stead, during non-conservative MT evolution, some frac¬ 
tion S of M p i,rlo is lost from the system and M* = 
— (1 — (5)M p i.rlo- Note that in none of the examples pre¬ 
sented here does the star fill its Roche-lobe during the 
orbital evolution calculation. 


3.4.1. Stable Conservative Mass Transfer 

During conservative MT, the evolution of J or b is gener¬ 
ally computed under the assumption that MT proceeds 
through an accretion disk. In the standard picture, the 
matter flowing through the inner Langrangian point ap¬ 
pears to be pushed into orbit about the host star by 
Coriolis forces (as viewed in the corotating frame). Sub¬ 
sequently, viscous str esses spread this material into a disk 
around the accretor dFrank ct al.lll985lh This disk trans¬ 
ports mass toward the accretor and angular momentum 
away from it. The latter is eventually returned to the 
orbit via torques operat ing between the donor a nd the 
outer edge of the disk (Lin fc Panaloizoul fT979h and a 
small residual angular momentum is transferred to the 
spin of the star. Our calculations with MESA include 
spin-up through accretion. However, as the material ac¬ 
creted by the star is only a tiny fraction of its total mass, 
spin-up due to accretion is negligible. 


3.4.2. Stable Non-Conservative Mass Transfer 

During n on-conservative MT, ch anges in J or b are com¬ 
puted as in lSoberman et all (11997T) 


Jorb \ / . \ 1 !‘l ATp]. RLO r- A/pi RLO 


(9) 


For comparison, the evolution of the orbital separation 
due to RLO is given by 


M 


RLO 


>o. (io) 

M p i 


This model assumes that a fraction 5 of M p i,rlo settles 
into a ring whose radius a r is a constant multiple y 2 of 
the orbital separation a. This mass is then lost from the 
system, taking with it the specific angular momentum 
of the ring. The remaining fraction of the mass (1 — 5) 
is assumed to be accreted onto the host star. Below we 
describe our choices for 7 and <5, based on the expected 
stability of MT. 


3.4.3. Considerations on the Stability of Mass Transfer 

The calculations presented in Section |4] show that con¬ 
servative MT is always stable. Instead, some care must 
be taken when considering non-conservative MT. 

For the systems under consideration, the stability of 
the MT phase depends mainly on the fraction of plan¬ 
etary mass leaving the system, its specific angular mo¬ 
mentum, and the response of the planet to mass loss. 
For a mass-losing planet the latter cannot be determined 
a priori and it is computed with MESA as the orbital 
evolution proceeds. We also have no prior knowledge 
of the values of 7 and 5 (the parameters regulating the 


amount of mass leaving the system and how much spe¬ 
cific angular momentum is carried away) that correspond 
to dynamical stability of MT. However, some guidance 
on the region of stability can be gained as follows. For 
systems with extreme mass ratios (q <C 1 , such as those 
considered here) and J 0 rb/Jorb given by Eqns. (0 and 
©■ the planetary mass during RLO changes according 
to (e.g.. iRanpaport et al.lll982r for a derivation) 


AI p i,rlo 

M p i 


( J orb | 

Jorb / tides 


A ( EeL\ _ -fcfpl,PE / 1 £ 

2 V^P 1 /therm M p‘ 6 2 


§-*y + § 

(ii) 


Here (R p i/R p i)therm is the fractional rate of change 
of the planet radius due to its thermal evolution and 
(Jorh/Jovh )tides is the fractional rate of change in or¬ 
bital angular momentum due to tides. Finally, £ = 
(dlni? p i/dlnAfpi) a d is the planet’s adiabatic logarithmic 
derivative of radius with respect to mass, not to be con¬ 
fused with the mass-radius relations in Figure [T| valid 
for thermal equilibrium models. A necessary condition 
for stable MT is that the denominator of Eqn. (fTTI) be 
positive. 

Our goal here i s to improve on the st able MT calcula¬ 
tions presented in lValsecchi et all (12014H by investigating 
the importance of irradiation effects and self-consistent 
models for the planet. To explore MT stability in the 
present work, we set 5 = 1 and consider 7 values for 
which the MT is expected to be dynamically stable. 
There is no loss of generality in this prescription since 
7 and 5 appear only as a product in Eqn. (1111) . We 
discuss possible realistic mass transfer configurations in 
Section0 while reserving a more detailed analysis of MT 
stability in hot-Jupiter systems to a future study. Such 
analysis should account for a broad range of 7 and 5 
values, as well as initial orbital configurations and prop¬ 
erties of the components (e.g., M*, M p 1 , and metallicity; 
see Section 0. 

With d = l and 7 ^ Owe are assuming that all RLO 
material leaves the system carrying away a specific an¬ 
gular momentum equal to 7 y/GM^a. The condition for 
stability, i.e., the requirement of a positive denominator 
in Eqn. m, then reduces to 



£ 

2 ‘ 


( 12 ) 


Eqn. 113 shows that an increase in the planetary radius 
with adiabatic mass loss (i.e., £ < 0 ) has a destabilizing 
effect, as expected. We test different values of 7 with 
MESA and find that, depending on the planetary core 
mass, M c , the computation becomes numerically difficult 
for 7 values higher than about 0.6 — 0.8. This suggests 
that the MT may indeed become dynamically unstable 
and that, according to the criterion in Eqn. m , £ is in¬ 
ferred to be close to zero. To avoid the instability region, 
while still considering somewhat significant angular mo¬ 
mentum loss from the system, we provide examples for 
non-conservative MT with 7 = 0.5 and 7 = 0.6 for the 
1 M® and 5 M® core cases, respectively, and 7 = 0.7 for 
the higher core masses. 
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4. RESULTS 

We consider a typical hot-Jupiter system comprising a 
1M© star and a lMj planet at solar metallicity. For the 
planet we took M c = 1 M®, 5 M®, 10 M®, 15 M® and, 
30 M®. The initial period was set to P 0 rb = 0.7 d in all 
cases. This corresponds to the orbital period at which a 
hot Jupiter with a 1.5-Rj radius would be at its Roche 
limit and it is chosen arbitrarily to have all systems start¬ 
ing with the same initial period “close enough” to the 
Roche limit. 

For the initial systems’ age we chose 2 Gyr (t — 
20 %tMs)> as it is at the low end of the ages of the 
currently known hot Jupiters closest to their Roche- 
limit separation (see Table 2 in lValsecchi fc Rasioll2014al 
and references therein; see also Section 0. Accordingly, 
we create 2 Gyr old stellar and planetary models with 
MESA’s “single-star” module. Each planet is irradiated 
according to the host star properties at 2 Gyr and the 
0.7d period (F* = 5.4 x 10 9 ergcm -2 s -1 ). 

The models are then used in MESA’s “binary” module 
to compute the orbital evolution. During this step both 
stellar and planetary evolution, as well as irradiation and 
photo-evaporation are computed self-consistently (Sec¬ 
tion 0 . For the stellar spin, we choose an initial value 
of fl* ~ 0.1 f 2 0 , where Q 0 is the orbital frequency. 
This is consistent with the observed slow stellar rota¬ 
tion rates for the tightest hot-Jupit er systems known 
(Q* ^ 0.1 — 0.2 0 o ; see Table 1 in IValsecchi fc Rasiol 
l2014alh Finally, as described in Section 0H we consider 
both conservative and non-conservative MT. 

The results are presented as follows. In Section 14.11 
we describe in detail the evolution of the 5 M® and 
30 M® planetary core cases. We take these as extreme 
examples among those considered here. In fact, as de¬ 
scribed below (Section 14. 2H , the behavior of the 1 M® 
core model needs a more in-depth investigation. For 
M c = 5 M® and M c = 30 M®, we mainly focus on the 
non-conservative MT evolution and briefly describe the 
conservative case at the end of each section. For these 
same core masses, we also investigate the effect of vary¬ 
ing the column density for irradiation in Section 14.1.31 
and the photo-evaporation recipe in Section 14.1.41 In 
Section 14.21 we present an overview of how the orbital 
period evolves as the planet loses mass, as well as how 
the planetary radius evolves with mass, for the full range 
of core masses. A summary of the results for all core 
masses and physical assumptions is presented in Table 0 
In all examples considered, the evolution is quite rapid 
after the planet is left with only a few percent of the 
envelope mass. Thus, we discuss the final stages of the 
planetary evolution qualitatively, guided by the relevant 
timescales entering the problem. 

4.1. Detailed Examples: Different Core Masses 
4.1.1. The 5 M® Core Model 

Figure 0 shows the evolution of a Jupiter with a 5 M® 
core undergoing non-conservative MT. The top panel dis¬ 
plays the evolution of various system and planetary prop¬ 
erties, while the bottom panel investigates a number of 
different timescales of the system. For convenience in 
studying the plot we arbitrarily subtract 2 Gyr from the 
time axis (this is just the time over which we evolved 
the star and the planet before inserting them into the 


binary evolution version of the MESA code; see discus¬ 
sion above). Hereafter, when describing various features 
in Figure 0 (as well as in Figure 0 we refer to the time 
marked on the axis, i.e., after subtracting off a 2 Gyr 
reference time. 

For the first ~70Myr the planet underfills its Roche 
lobe while the orbit shrinks due to the tides transferring 
angular momentum from the orbit to the stellar spin, in 
a vain attempt to try to spin-up the star. Eventually, 
when P or b shrinks to ~ 0.5 d (a ~ 0.01 AU) the planet 
fills its Roche-lobe. This occurrence can be seen when the 
dotted curve in Fig. 0 (indicating the Roche lobe radius) 
merges with the planetary radius curve (solid blue). 

The system remains in Roche lobe contact for the en¬ 
suing 1.6 Gyr, while the planet loses more than 95% of 
its mass. Concomitantly, the orbital period grows from 
about 0.5 to 3.9 days. During this same interval, the ef¬ 
fective temperature of the planet decreases from about 
2500 K to about 1300 K, due largely to the increasing 
distance between the host star and the planet. 

During the RLO phase there is a competition between 
the tidal effects which tend to drive orbital decay and 
mass transfer which generally drives orbital expansion 
[see Eqn. Col)]. Both effects operate simultaneously, and 
the one that dominates depends on the response of the 
planet’s radius to mass loss. For Roche-lobe overflowing 
objects with mass lower than the accreting star, there is 
a universal relation between the orbital peri od, and mass 
and r adius of the donor [see, e.g., Eqn. (2) of lHowcll ct all 

l200i| : 

P mh ~ 0.4 (M vl /Mj)- l /\R pX /Rjf/ 2 days (13) 

Thus, if the mass transfer remains stable, then the or¬ 
bital period will grow or shrink in accordance with how 
the planetary radius changes with mass loss. In this ex¬ 
ample with a 5 M® core, we can see from Fig.0that the 
planet’s radius does not vary very much during most of 
the RLO phase, and in fact, slightly increases to 1.6 Rj. 

Therefore, P or b will simply grow mostly as oc M pl ' . By 
the end of the RLO phase the planet’s mass has shrunk 
to 8.5 M®, and, according to Eqn. m P or b should be 
about 4 days, which it is. The Roche lobe overflow phase 
ends when the combination of tidal decay of the orbit 
and planetary expansion is no longer sufficient to keep 
the planet filling its Roche lobe, and the planet, which 
is now quite low in mass, starts to shrink well inside its 
Roche lobe. 

The post-RLO phase is driven mainly by photo¬ 
evaporation. In fact, the tidal decay timescale becomes 
longer than a Hubble time for the majority of this phase 
due to the longer orbital period and the greatly reduced 
planetary mass. Photo-evaporation removes mass from 
the planet at a rate of about 10 -16 — 10 _ 15 M©yr _1 
( 10 10 — 10 11 g s -1 ) and it does not affect the orbital sep¬ 
aration significantly [see Eqn. 0]- For the 5.4 Gyr dura¬ 
tion of the post-RLO phase the planet remains in a 3.9 d 
orbit and its mass decreases from about 8.5 M® to about 
5.4 M® (when M env /M p \ ~7%), when the calculation 
ends because of convergence problems. 

The bottom panel of Fig. 0 displays various timescales 
for this evolution, including the mass loss timescales via 
the photo-evaporative wind, |M p i/Mpi|pE, and due to 
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t-2(Gyr) 

Figure 3. Detailed evolution of some of the components and orbital parameters (top) and of the relevant timescales (bottom) for a Jupiter 
with M c = 5M®. For this evolution we adopt the parameters 5 = 1 and 7 = 0.6. The subscripts “RLO”, “PE”, and “Tides” refer to 
the timescales associated with mass loss due to Roche-lobe overflow and photo-evaporation, and to tidal decay, respectively. For the latter 
we note that only stellar tides play a significant role, as the planetary tides timescale is always longer than a Hubble time. The peaks in 
the tidal timescales occur when the contributions from stellar and planetary tides cancel out. The dotted blue curve in the upper panel 
indicates the Roche-lobe radius, while the dotted curve in the bottom panel is the stellar nuclear timescale. The RLO phase ends after 
about 1.6 Gyr, when the system is about 3.6 Gyr old. For clarity, the timescale for the evolution of the planetary radius (bottom panel, 
green line) is computed taking the median of 20 consecutive values. The calculation ends when M en v/M p 1 < 7% because of convergence 
problems. 


RLO, |Afpi/M p i|RLO, as well as the timescale for tidal 
decay of the orbit, | J 0 rb/Jorb \tides , and for the thermal 
expansion/contraction of the planetary radius, |i? p i/i? p i|. 
What we see is that, after RLO commences, the black 
curve, showing the mass transfer timescale (due to RLO), 
lies a factor of about 8 lower than the tidal decay 
timescale, and much lower than the mass loss timescale 
associated with PE winds. We learn from Eqn. (fTTl) that 
the mass loss rate due to RLO is therefore essentially pro¬ 
portional to the rate of decay of the orbit due to tides, 
and inversely proportional to the denominator which is 
£/2 — 7<5 + 5/6 = £/2 + 0.23. Putting these together 
implies that the denominator must equal approximately 
1/8. From this we can infer that the effective adiabatic 
index of the planet during most of the RLO phase is 
£ ~ —0.2, in other words the planet reacts to adiabatic 
mass loss by slightly expanding. Later in the RLO phase, 
the tidal decay timescale greatly increases due to the in¬ 
creasing orbital separation, but the thermal expansion 
timescale of the planet decreases dramatically to pick 


up the slack of the declining tidal effects. We also see 
that the mass loss rate in a PE driven wind dramati¬ 
cally increases (i.e., the timescale decreases) due to the 
increasing radius and the decreasing mass of the planet 
[see Eqn. (©]. 

Even though our evolution calculation ends when the 
envelope mass fraction drops to about 7%, we argue that 
eventually, the host star will approach the end of its main 
sequence lifetime and it will begin expanding. The in¬ 
crease in A* will cause a concomitant increase in both 
the tidal decay rate and in the amount of irradiation 
received by the planet. The latter occurs because A* 
evolves faster than T* in Eqn. ©. Such an increase in 
irradiation will yield an increase in A p i and, as a conse¬ 
quence, even stronger photo-evaporation. What happens 
to this system after the envelope of the planet has been 
removed, leaving only its rocky core, depends mainly on 
the tidal timescale compared with the stellar nuclear evo¬ 
lution timescale. The former determines the rate of tidal 
decay (leading to a second planetary RLO), while the lat¬ 
ter determines the rate of stellar expansion (leading to 
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stellar RLO). Given the (very steep) dependence of the 
tidal d ecay timescale as (a/i ?*) 8 fe-g.. I Valsecchi fe Rasiol 
1201414) . we argue that tides will likely drive the planet 
down to its Roche limit. Here, it will be stripped of 
the remaining envelope and, assuming a constant den¬ 
sity for the core, it will be consumed at approximately 
constant P nr b on the ra pi d tida l decay timescale. This 
follows from iPaczvriskiD s ( 1971 1 approximation for the 
Roche limit separation or = (i? p i/0.462)(M*/M p i) 1 / 3 oc 
(Mt/ppi) 1 ' 3 , where p p 1 denotes the mean density of the 
planet. 



Figure 4. Same as Figure [ 3 ] but for a Jupiter with M c = 30 M® 
and 7 = 0.7. For clarity, |M p i/M p i|rlo is computed by taking 
the median of 10 consecutive values. 


The above discussion was for the case of non¬ 
conservative MT (with 7 = 0.6). In the case of conser¬ 
vative mass transfer, we find that the orbital evolution 
proceeds similarly to the examples presented above, but 
the duration of the various phases is different. In fact, 
when MT is conservative, the tidally-driven orbital decay 
is counteracted by a more substantial RLO-driven orbital 
expansion [<5 = 0 in Eqn. (1101) ]. which results in a slower 
overall evolution. This is summarized in Table [5] In par¬ 
ticular, the RLO phase lasts for about 4.5 Gyr, leaving 
an 8.3 M® planet in a 3.4 d orbit. The remainder of the 
evolution proceeds similarly to the non-conservative case. 
However, the evolution when little envelope mass is left 
is shorter. This is due to the irradiation-driven decrease 


in planetary mass loss timescale (|M p i/M p ppe|) due to 
the star approaching the end of its main sequence. 

4.1.2. The 30 JW® Core Model 

Figure Q] shows the evolution of a Jupiter with a 30 M® 
core undergoing non-conservative MT (with 7 = 0.7). As 
in Fig. [3] (for the case of a 5 M® core), the top panel 
presents the evolution of various system and planetary 
properties, while the bottom panel displays a number of 
different timescales of the system. For the first ~ 70Myr 
tides cause the orbit to shrink and the planet fills its 
Roche lobe when P or b — 0.5 d (the same as for the 5 
M® core case). 

At the onset of RLO, the orbit expands as mass 
is removed from the planet at a rate of ~ 10 -13 — 
10 ~ 12 Mq yr _1 (10 13 — 10 14 g s _1 ). However, after about 
a Gyr, the orbit begins shrinking once the period has 
grown to only ~0.8 days. Overall, we follow the RLO 
phase for about 2.1 Gyr, to the point where the plan¬ 
etary envelope has been completely removed. At this 
time, the orbital period has shrunk to 0.3 d, and the ef¬ 
fective temperature of the planet has increased to nearly 
3000 K. 

The difference in evolutionary history between this 
planet with a 30 M® core and the one with a 5 M® core 
(described in Section 14.1.11) results largely from the fact 
that planets with more massive cores exhibit an earlier 
(i.e., at higher M p 1 ; Figure [HD decrease in planetary ra¬ 
dius with continuing mass loss. During the RLO phase 
with a 30 M® core, the mass of the planet decays essen¬ 
tially as a power law in time, while the radius does not 
decrease significantly until M p \ drops below ~60 M®. 
From Eqn. (fl3l) we can deduce that this combination of 
mass and radius changes will lead to a steady increase 
in the orbital period. However, once the radius starts to 
decline substantially, the R p ^ dependence in Eqn. m 
dominates the orbital period evolution, and P 0 rb starts 
to decay. 

We can gain some further insight into the evolu¬ 
tion of this system by considering the timescales dis¬ 
played in the bottom panel of Figure |4] For most of 
the evolution, the tidal decay timescale (| J or b/Jorb| tides) 
is about 6 times longer than the mass loss timescale 
(|M p i/M p i|rlo)- Furthermore, at least for the earlier 
portion of the RLO phase, both the contraction timescale 
of the planet (|i? p i/i? p i|) and the timescale of photo¬ 
evaporation (|M p i/M p i|pe) are both longer yet than the 
tidal decay timescale. From this we can infer that the 
denominator in Eqn. ED must be ~ 1/6. In turn, we 
can conclude that £ ac j ~ 0.07 (i.e., very close to zero). 
During this earlier portion of the RLO evolution, the or¬ 
bital period grows and the mass of the planet declines, a 
combination that leads to an ever increasing tidal decay 
timescale (i.e., weakening of the tidal evolution of the 
orbit). 

During the later portion of the RLO phase, both 
|-Rpi/-Rpi| and |M p i/M p i|pe become shorter than the tidal 
driving timescale. This means that there is close com¬ 
petition in the numerator of Eqn. ED for maintaining 
Roche-lobe contact among all three terms: (i) tidal de¬ 
cay; (ii) the photo-evaporative mass loss term (tending 
to shrink the Roche-lobe radius of the planet and main- 
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tain RLO), and (iii) the shrinkage of the planet as it 
loses mass, which tends to push the planet back within 
its Roche-lobe. Apparently, the combination is sufficient 
to maintain RLO until we terminate the evolution. 

Our calculation ends when the planetary envelope has 
been completely removed. However, the core itself may 
well eventually undergo RLO and, if the core has a con¬ 
stant density, it will be consumed at constant orbital 
period on the rapid tidal decay timescale (< 100 Myr, at 
the end of the evolution in Figure |4j). 

The conservative mass transfer case for the same model 
with a 30 M® core differs only insofar as the duration 
of the RLO phase is concerned (see Table © ■ This is 
the same as what we found for the conservative vs. non¬ 
conservative cases with a 5 M® core. In fact, during con¬ 
servative MT, the RLO phase until the envelope is re¬ 
moved lasts about 3.8 Gyr. At the end of the calculation 
Porb = 0.3 d. 

4.1.3. Varying the Column Density for Irradiation 

For highly irradiated planets with M n i < Mj, the 
q uasi- equilibrium mass-radius relations of lFortnev et al.1 
(l2007ll in Figure Q] are bracketed by those computed 
with MESA for irradiation absorption column densities 
of S p i = (1—100) g/cm 2 . We tested the effect on plane¬ 
tary RLO evolution by increasing E p i by two orders of 
magnitude in the 5 M® and 30 M® core mass models. 
We find that the overall evolution does not change sig¬ 
nificantly (see Table ©. Qualitatively, for the 5 M® core 
case, at the beginning of the calculation the radius at 
higher E p i is a few percent larger. As a result, RLO starts 
at a longer orbital period and, for the same evolutionary 
time, it continues at longer P or b- When the planet re¬ 
treats back within its Roche-lobe, the longer orbital pe¬ 
riod yields a less severe mass loss via photo-evaporation. 
Consequently, the planet retains a higher envelope mass 
fraction (/ in Tabic [2]) for a longer time. A similar be¬ 
havior occurs for the 30 M® core model, but the percent 
difference in radius (orbital period) increases from ~ 5% 
to ~ 15% (from ~ 5% to ~ 20%) between the begin¬ 
ning and the end of the calculation. Again, mass loss 
proceeds on a longer timescale and, for the same evolu¬ 
tionary time, the models with higher E p i retain a larger 
fraction of the envelope for longer. 

4.1.4. Varying Photo-Evaporation 

Formally, the prescription of lMurrav-Clav et all (j2009l ) 
uses Pxuv ~ 10 4 erg cm -2 s -1 as a threshold between 
the energy limited and radiation/recombination limited 
regimes. Below we investigate whether our results for the 
5 M® and 30 M® core cases change significantly by using 
Eqn. © when Pxuv > 10 4 erg cm -2 s _1 and Eqn. 0 
otherwise. The results are summarized in Table [2] (de¬ 
noted with “e-Lirn + rr-Lim” in the column named 
“PE”). 

The planetary mass loss rate in the radia¬ 
tion/recombination limited regime is slower. This 
naturally leads to a longer duration of the RLO phase 
for all examples considered. For the 5 M® core under¬ 
going conservative MT, R p \ is a few percent larger. As 
a result, the orbit evolves at longer period. However, 
because of the longer RLO phase, by the time the planet 
is left with about 10% of its envelope mass the star is 


approaching the end of its main sequence. Consequently, 
the increase in R p \ driven by the increase in stellar irra¬ 
diation leads to faster mass loss via photo-evaporation 
towards the end of the calculation. For the 5M® core 
undergoing non-conservative MT, the evolution follows 
this same line of logic, but the planetary radius is a few 
percent smaller for most of the evolution. As a result, 
the orbit evolves at shorter period. For the 30M® core 
case, the response of R p \ to mass loss for the different 
photo-evaporation prescriptions differs by only 1 — 2%. 
Therefore, only the duration of the RLO phase changes 
significantly, while the orbital period at the various 
stages of Table [2] is not significantly affected. 

In summary, none of our basic results would change 
significantly if we were to utilize both the energy-limited 
and radiation/recombination-limited prescriptions for 
photo-evaporation over the entire evolution. However, 
certainly a number of the details of the evolutionary mod¬ 
els would look somewhat different. 

4.2. A Range of Evolutionary Models 

A range of evolutionary models covering different core 
masses for both conservative and non-conservative MT 
are presented here. The results are summarized in Fig¬ 
ures [5] and [© respectively. The top panels show the evo¬ 
lution of the orbital period with the planet’s mass. For 
comparison, the positions of the observed planets are the 
orange open circles. The bottom panels show the corre¬ 
sponding evolutionary tracks for planetary mass and ra¬ 
dius. The evolution tracks for different core masses are 
denoted with different colors and line-styles. For all cases 
considered, the MT proceeds on a timescale longer than 
the thermal timescale of the planet. Thus, the planet re¬ 
mains in near thermal equilibrium throughout the MT. 
Here we note that our 1M® core model shows a severe 
increase in R p 1 that needs further investigation. In fact, 
R p i increases up to about 10 Pj when M p \ ~ 1M®. 
However , we note also t hat t he core-less models in Fig¬ 
ure 8 of lFortnev et all (120071 ) have radii ranging up to 
P p i ~ 2.3 Pi w hen M n 1 ~ 30M® (the lowest mass con¬ 
sidered by iFortnev et all (20071 for a core-less Jupiter), 
depending on the age of the planet and the level of ir¬ 
radiation. This is consistent with our radii for masses 
down to ~10 M®. We omit the 1M® core model from 
the subsequent discussion. 

The evolution in time in the M p \ — P or b plane is down¬ 
ward (i.e., shrinking mass), initially toward longer peri¬ 
ods, but then decaying toward shorter P or b- The tracks 
in the (P p i — M p \) diagram proceed from the right to 
the left. The shape of the evolutionary tracks in such 
diagrams depends on the competing effects of stellar 
tides (tending to shrink the orbit) and RLO (tending 
to expand the orbit). As described in Section 14.1.11 the 
P p i — M p i tracks in the bottom panels of Figures [5] and 
[G] partly determine which mechanism ends up dominat¬ 
ing in terms of net orbital contraction or expansion [see 
Eqn. mi 

At the onset of RLO, the nearly constant or increasing 
P p i with decreasing M p \ causes the RLO term to dom¬ 
inate over the tidal contribution. As a result the orbit 
expands. This behavior persists until the mass-radius 
relation begins steepening, becoming more positive. At 
this point, the tidal term becomes more significant than 
the RLO term and the orbit begins to shrink. In some 
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Figure 5. Planetary mass as a function of the orbital period 
(top) and mass-radius relation (bottom) for conservative MT (with 
5 = 0). Because some of the evolutionary calculations become nu¬ 
merically challenging when the mass fraction in the envelope drops 
below a few percent, we show the evolution up to when M en ^/M p \ 
drops just below 5%. The open orange circles are confirmed ex¬ 
oplanetary systems (NASA Exoplanet Archive, 13 January 2015) 
with observationally inferred M p \ and P Q rb » and hosting one planet 
only, for simplicity. The system GJ 436b dButler et al 1200411 would 
be located at P 0 rb — 2.6 d and M p \ ~ 22 M®. The colored lines are 
our evolutionary models for different core masses. The colored ar¬ 
rows along each evolutionary track in the top panel mark 1 Gyr in¬ 
tervals and denote time evolution. For the remaining symbols: “xs” 
mark the end of RLO, “ds” mark times when the mass in the en¬ 
velope drops below / = 30%, and 5% of the total mass. For the 
5 M® core model, the various symbols in the top panel are all su¬ 
perimposed at the end of the evolution, where the planet spends a 
few Gyrs. 

cases (M c = 5Mg, 10M@, 15M®), the combination of 
the planet shrinking in response to mass loss and the 
weakening tidal forces with decreasing M p \ causes the 
RLO phase to terminate and the planet to detach (“x” 
symbols). In a few cases (M c = 5M®,10M® during 
conservative MT), as the star approaches the end of its 
main sequence, the increasing amount of irradiation re¬ 
ceived by the planet causes R n ] to increase again (e.g., 
iFortnev fe Ncttelmannll2010f ). Note that none of the ob¬ 
served exoplanets in Figures [5] and M is currently in 
RLO. Therefore, observations should be compared with 
the portion of the evolutionary tracks where the planet 
is detached. There are some 6-7 observed systems shown 
in these figures which are in the vicinity of our evolu¬ 



Figure 6 . Same as Figure l5l but for non-conservative MT (with 
5=1). We use 7 =0.5, 0.6, and 0.7 for M c = 1 Af®, 5M®, 
and > 10 M®, respectively. Note, for the 5M® core model the 
evolution terminates because of convergence problems when the 
envelope contains ~ 6 % of the total mass. For the 5 M® core 
model described in detail in the text, the various symbols in the 
top panel are all superimposed at the end of the evolution, where 
the planet spends few Gyrs. 

tion tracks, after Roche-lobe overflow has stopped, and 
our mod els may be direct l y app licable to them. Note 
also that lEhrenreich et al.1 (l2015h recently reported the 
discovery of a large exospheric clou d surround i ng th e 
Neptune-mass exoplanet GJ 436 ( Butler et al.1 12004 1. 
composed mainly of hydrogen atoms. The average ob¬ 
served mass loss rate implies an efficiency for converting 
X-ray and extreme UV energy into mass loss of about 
1%. In Figures [5] and [6] this planet would be located at 
P or b — 2.6 d and M p \ ~22M®. 

Finally, we note that the shape of the evolutionary 
tracks does not change significantly in the M p \ — P or b di¬ 
agram between conservative and non-conservative MT. 
However, as summarized in Table [2] and explained in 
Section ITU whether mass is lost from the system or not 
does affect the duration of the various phases mentioned 
above. 

5. DISCUSSION 

Our calculations seem very promising for explaining 
some of the super-Earth and sub-Neptunes-type plan¬ 
ets whose bulk density suggests that they consist of 


Rpi(Rj) 
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a core (rocky or icy) surrounded by a H/He envelope 
comprising; up to tens of a percent of the total mass 
(jLopez &; Fortnevll2014lh While this agreement with our 
MT models is very encouraging, our calculation neglects 
some important effects which we discuss below. Specif¬ 
ically, in Sectio n Hull we discuss possible MT scenarios, 
while in Section 15721 we discuss our assumptions on plan¬ 
etary tide s. A discussion of observational signatures is 
in Section IQl 

5.1. Mass Transfer Scenarios and Stability 

In this work we neglected the effects of magnetic fields 
and stellar winds on the RLO material. Th ese mecha¬ 
nism s can affect the f l ow of MT (c.g., iCohen fc Gloced 
l 2012 t I Owen k, Adams! I2014T) , potentially playing a cru¬ 
cial role in determining whether any mass is transferred 
or an accretion disk ever forms. On the opposite side of 
conservative MT, is the case where mass is blown away 
directly from the planet or the inner Lagrangian point 
(Li). For a lMj plane t and a 1 M& st ar, L\ is located at 
~ 0.93 a [Eqn. (10) of lLai ct alJl20ldj |. In the formalism 
of Section [3] these scenarios can be reproduced by setting 
S = 1 and 7 = 1 (mass lost from the planet) or 7 = 0.97 
(mass loss from Lf). In this configuration, Eqn. m re¬ 
quires £ to be larger than (0.27-0.33) in order for the 
MT to be dynamically stable. The values of £ inferred 
in Sections 14.1.11 and 14.1.21 for our extreme core masses 
(£ = —0.2 and 0.07 for the 5M® and 30M®, respectively) 
suggest that the MT will likely be dynamically unstable 
if 7 > 0.97. Furthermore, we have performed test runs 
with MESA fixing S = 1 while varying 7 . We find that 
the MT may become dynamically unstable for 7 > 0.8. 
In fact, the computation becomes numerically difficult, 
with the integrator time-step dropping to less than one 
month. This limit on 7 for stability suggests that the ef¬ 
fective value of £ is close to zero, or negative (as directly 
computed within MESA). 

To gain better insight into more physically meaning¬ 
ful values of the parameter 7 , consider the case where 
matter flows in a narrow stream from the L\ point, and 
forms a ring around the star which is then blown away. 
This case has been discussed in the c ontext of MT sta¬ 
bility in stellar binary systems (e .g., iHut fc Paczvnskil 
11984 Verbunt fc Rappaportlll98§ h If, in fact, the life¬ 
time of the gas in the ring, before being blown away, 
is shorter than the viscous timescale for the ring, then 
the angular momentum may not be returned to the or¬ 
bit. In this case, the angular momentum leaving the 
system would be determined by how much angular mo¬ 
mentum a particle has before it is blown away (i.e., by the 
mean ring rad i us rn ) . We combine calculations of by 
iLubow fc Shul (119751) and lHut fc Paczvnskil (|1984f) . cov¬ 
ering Mpi/iW* down to 10 —3 , with our own calculations, 
extending these down to M v \/M t = 10” 6 . We find r d 
to range between ~ 0.7 — 0.9 for a 1 Mj — 10 M® planet. 
This range implies 7 values between ~ 0.84 — 0.95 and 
a value of £ for stability between ~0—0.23 [Eqn. (fl2l) ]. 
This scenario appears borderline between stable and un¬ 
stable MT. Clearly any scenario where 7 > 1 would be 
dynamically unstable. 

If the MT were indeed dynamically unstable, one could 
envision the system undergoing a common-envelope-like 
evolution (‘CE’), that is somewhat distinct from the stan¬ 


dard binary stellar evolution picture (IWebbinkl Il984f) . 
We envision that the envelope of the planet would quickly 
flow though the inner Lagrange point and form a disk-like 
structure around the host star. The core of the planet 
would then find itself orbiting within this ring or disk 
of envelope material. Depending on the detailed core¬ 
disk interactions, the core may spiral-in toward the host 
star and eject the disk material. This is distinct from 
the usual CE scenario in that the envelope of the planet 
becomes bound to the host star, rather than the planet, 
and it is the planet which ejects its own remnant enve¬ 
lope via tidal and viscous interactions. Insights into the 
outcome of such a phase can then be gained from the 
energy equation 


O-CE 


GM*M C 

2af 


GM*M C 

2 ai 


GM,tf mv 

2 a; 


(14) 


Here, a®E represents the efficiency with which the planet 
core’s orbital energy can be used to unbind the envelope 
material which is now in a ring around the host star. The 
parameters a, and af denote the orbital separation at 
the onset and at the end of this ‘CE’ phase, respectively. 
Eqn. m can be directly solved for the ratio of a//a*, 
and we find: 

aj_ _ «CeM c 

cij a ceM c + M env 

For plausible values of a ®e near unity, this expression 
yields af/a,i ~ M c /M p \. This, in turn, implies that the 
orbital separation would decay by a factor of more than 
an order of magnitude if the mass transfer is unstable at 
the onset of RLO. 

Thus, all cores considered in this work would reach 
their own Roche limit in the event of unstable RLC0. 
A new RLO phase would then begin driven, this time, 
by the interplay of viscous drag on the core due to the 
remaining envelope mass orbiting the host, and the back- 
reaction from continuing (now stable) MT through RLO. 
The former would tend to shrink the orbit, while the 
latter would cause the orbit to expand. 

The actual outcome of an unstable MT phase from a 
hot Jupiter to its stellar host is still an unexplored ques¬ 
tion which we are currently addressing via smoothed- 
particle hydrodynamics simulations. Furthermore, mag¬ 
netohydrodynamics (MHD) simulations are needed to 
determine the role played by ma gnetic fields. Here we 
note that iMatsakos et al.l (|2015l) recently investigated 
magnetized star-planet interactions via 3D MHD nu¬ 
merical simulations including stellar wind and photo- 
evaporative mass loss. They found that, depending on 
the planet’s magnetic field and outflow rate, as well as 
the stellar gravitational field, the star can accrete part 
of the mass lost by the planet via photo-evaporation. 
iPillitteri et al.1 (1201511 proved this to be a plausible sce¬ 
nario via far-ultraviolet observations of HD 189733. 


5.2. Planetary Tides 

In this work we have taken planetary tides to be effi¬ 
cient in keeping the spin of the planet tidally locked. This 


8 If it* is too large, the planet-core will plunge into its atmo- 
sphere before filling its Roche lobe; the critical p er iod for tida l 
breakup of a rocky body can be <5 hours ijRappaport et aJ1l2Q13l ). 
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assumption is justified by the magnitude of the tidal syn¬ 
chronization timescales discussed in Section [3] However, 
we neglected the resulting effects of tidal dissipation on 
the planetary structure. These include potentially sig¬ 
nificant heating, inflation, and resultant stronger mass 
loss. To gain some intuition into their significance we 
compute the power deposited between two consecutive 
integration steps and compare it with the planet lumi¬ 
nosity, L p i, at each integration step. For the former we 
use Ltide = 5 / p i|flo(t + A t) — flo(f)|/A i, where / p i is the 
moment of inertia. We find that L t ide/Api decreases from 
1CT 4 to 10 -11 throughout the calculation. This suggests 
that the power deposited in the planet via tides may not 
affect its structure dramatically. 

5.3. Observational Signatures 

As discussed in IValsecchi et al.l (12014 1. this evolu¬ 
tionary scenario has several observational consequences. 
First, the properties of stars hosting hot Jupiters and 
those hosting Super-Earth and mini-Neptune-type plan¬ 
ets should be similar. For example, if all RLO material 
is lost from the system, the host stars should have sim¬ 
ilar mass and composition, but super-Earth and mini- 
Neptune host stars should be somewhat more evolved. 
Furthermore, depending on the core mass and the de¬ 
tails of the MT process (whether any mass is lost from 
the system and the location where angular momentum 
is removed), a system may spend enough time in RLO 
that it might be possible to observe planets in such a 
phase. If MT really proceeds through an accretion disk, 
this may produce observational signatures (e.g., line ab¬ 
sorption of stellar radi ation and time- dependent obscu¬ 
ration of the starlight; iLa.i et al.ll20T7)lb Finally, the re¬ 
sults in Figures [5] and El and in Table [2] suggest that, if 
this model is a viable formation channel for super-Earth 
and mini-Neptune-type planets, there should be a cor¬ 
relation between M p \ and P 0 rb- Specifically, the more 
massive planets should be found at shorter orbital peri¬ 
ods. This is a major result when c ompared to the sim¬ 
ple model of IValsecchi et all ( 2014 1. In fact, without a 
sclf-consistent calculation of planetary evolution and ir¬ 
radiation effects for a spectrum of core masses, we had 
found no trend between final M p \ — P n rh pairs. T he d if¬ 
ferent selection of core masses adopted in IValsecchi et al.l 
( 2014 1 does not allow for a one-to-one comparison of 
those results with the findings of this paper. However 
these different results can be understood simply via the 
way in which the orbital separation evolves in response 
to planetary mass loss, and via the incomplete descrip- 
tio n of the pre-det ermined mass-radius relations utilized 
in IValsecchi et al.l (|2014f l. The latter fixed M p \ — R p \ re¬ 
lations are shown in Figure [7] compared with those com¬ 
puted with MESA during the course of the evolutions. 

This inverse correlation between orbital period and 
planetary mass that we have found also suggests that 
dynamically stable MT phases like those presented 
here do not seem to represent a viable channel for the 
formation of the so-called bil tra-short-p eriod planets’ 
(e.g., ‘USPs’; iSanchis-Oieda et all 120141 and reference 
therein). These planets have typical radii smaller 
than 2 R m (correspon ding to a mass of about 5 M®; 

I Weiss &; Marcvl[2014f) and orbital periods shorter than 



Figure 7. Mass-radius relations. In black are the planetary 
models computed in this work, while in red are those used in 
IValsecchi et al.l (|2014l b 

one day. 


6 . CONCLUSIONS 

In this paper we have presented the first evolution¬ 
ary calculations of irradiated hot-Jupiters undergoing 
tidally-driven Roche-lobe overflow (RLO). We found 
that, depending on the size of the planetary core and 
the details of the mass transfer, the RLO phase and, in 
turn, the detached phase (after RLO had ceased) can 
last from a few to several Gyrs. For the smaller core 
masses (M c < 15 M®), after most of the envelope has 
been removed during RLO, the planet spends a few Gyrs 
losing mass via photo-evaporation at nearly constant or¬ 
bital period. This is consistent with the density of known 
supcr -Earths and sub-Neptun es, for which detailed mod¬ 
eling (ILqdcz fc Fortne-(l2014f ) place tens of percent of the 
total mass in a H/He envelope surrounding a rocky core. 

As noted above, we find an inverse correlation between 
the core mass of the planet and its final orbital period. 
This results from the basic fact that, in general, the ir¬ 
radiated planets with larger core masses decrease in ra¬ 
dius faster/earlier with mass loss than planets with lower- 
mass cores. Final orbital periods of < 1 day appear to 
require large core masses of > 15 M®, which are likely 
substantially hig her than the masse s of the USP planets 
found bv lSanchis-Oieda et all (120141) . Thus, the scenario 
we are presenting here probably does not account for the 
USPs. 

In this work we have considered a coarse grid of core 
masses and one initial binary config uration . H owever, 
as sum marized in Tables 1 and 2 of IValsecchi fe Rasiol 
l2014al the tightest observed hot-Jupiter systems com¬ 
prise a variety of stellar and planetary masses (~ 0.87 — 
1.33 Mq and ~ 0.46 — 1.49Mj), as well as metallicities 
(. Fe/H ~ —0.35 — 0.22) and ages (1.5—13Gyr). All these 
parameters may affect the efficiency of tides and, thus, 
the evolution of a Jupiter undergoing RLO. This vari¬ 
ety of properties requires exploration of a more refined 
grid in parameter space of initial component and orbital 
properties, as well as planetary core masses (especially 
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in the M c < 15 M 0 regime; Figures [5] and [6]). Such an 
extended parameter space study should also explore in 
more detail the boundaries for MT stability. Future ob¬ 
servations of increasingly massive super-Earths and sub- 
Neptune-type planets in increasingly tighter orbits might 
provide an important observational test of the ideas pre¬ 
sented here. 

Finally, we remark that the value of jS, the parameter 
describing mass and specific angular momentum loss, is 
crucial to the stability of RLO mass transfer. Where 
matter goes after RLO, and how much of it is actually 
accreted by the host star or ejected from the system, can 
only be determined by hydrodynamic calculations. The 
results of such calculations could well determine whether 
hot Jupiters can undergo the kind of stable RLO mass 
transfer described in this work. 
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Table 2 

Summary of Results. 


M c 

7 

S p i 

PE 

AtRLO 

P Q rb at the 

/ at the 

At when 

P orb when 

At when 

P n rh when 

t when 






end of RLO 

end of RLO 

/ ~ 20% to 30% 

/ ~20% 

/ ~7% to 20% 

/ ~7% 

/ ~7% 

(M®) 


(gem -2 ) 


(Gyr) 

(d) 

(%) 

(Gyr) 

(d) 

(Gyr) 

(d) 

(*MS) 

5 

- 

1 

e-lim 

4.5 

3.4 

39.8 

0.2 

3.4 

3.1 

3.4 

1.0 

10 

- 

1 

e-lim 

4.4 

1.6 

43.1 

0.5 

1.6 

2.9 

1.6 

1.0 

15 

- 

1 

e-lim 

4.2 

1.2 

43.1 

0.6 

1.1 

2.0 

0.9 

0.9 

30 

- 

1 

e-lim 

3.7 

0.5 

7.0 

0.3 

0.6 

0.3 

0.5 

0.6 

5 

0.6 

1 

e-lim 

1.6 

3.9 

41.1 

0.1 

3.9 

5.2 

3.9 

0.9 

10 

0.7 

1 

e-lim 

1.2 

1.7 

46.0 

0.3 

1.7 

3.4 

1.7 

0.7 

15 

0.7 

1 

e-lim 

1.3 

1.2 

48.2 

0.5 

1.2 

3.4 

1.1 

0.8 

30 

0.7 

1 

e-lim 

2.0 

0.5 

7.0 

0.3 

0.6 

0.4 

0.5 

0.4 

5 

- 

1 

e-lim + rr-lim 

6.9 

3.8 

37.7 

0.2 

3.8 

1.6 

3.8 

1.1 

30 

- 

1 

e-lim + rr-lim 

5.6 

0.5 

7.0 

0.4 

0.6 

0.3 

0.5 

0.8 

5 

0.6 

1 

e-lim + rr-lim 

3.7 

3.5 

40.1 

0.2 

3.5 

3.8 

3.5 

1.0 

30 

0.7 

1 

e-lim + rr-lim 

3.0 

0.5 

7.0 

0.4 

0.6 

0.3 

0.5 

0.5 

5 

- 

100 

e-lim 

4.7 

3.7 

40.2 

0.2 

3.7 

3.5 

3.7 

1.0 

30 

- 

100 

e-lim 

4.2 

0.5 

7.0 

0.4 

0.7 

0.4 

0.5 

0.6 

5 

0.6 

100 

e-lim 

1.7 

4.3 

41.7 

0.2 

4.3 

5.3 

4.3 

1.0 

30 

0.7 

100 

e-lim 

2.0 

0.7 

19.4 

0.3 

0.7 

0.6 

0.5 

0.5 


Note. — The reshlts include / = M enw /M p \ values down to 7% (but for the values in boldface where / >8%), to allow a comparison between different MT assumptions and M c 
values. In fact, for some of the models considered MESA does not converge when M env /M p \ drops below about 7%. Time intervals are denoted with At. Column PE denotes the 
photo-evaporation prescription, with e — lim and rr — lim indicating the energy-limited and radiation/recombination-limited regimes, respectively. The parameter t denote the age of 
the system in units of the stellar main sequence lifetime. The first four examples are for conservative MT. For non-conservative MT we use <5=1 and vary 7 as summarized in the 
second column. Note, the 30 M® core cases never detaches, but for the last case (non-conservative MT evolution with E p il00g cm -2 ). In this case, 0.5 Gyr after the first 2 Gyr long 
RLO phase, the planet fills its Roche lobe again). 
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